Dissecting the Cygnus Region with TeV Gamma Rays and Neutrinos 
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Recent Milagro observations of the Cygnus region have revealed both diffuse TeV gamma-ray emission and 
a bright and extended TeV source, MGRO J2019+37, which seems to lack an obvious counterpart at other 
wavelengths. Additional study of this curious object also promises to provide important clues concerning one of 
the Milky Way's most active environments. We point out some of the principal facts involved by following three 
modes of attack. First, to gain insight into this mysterious source, we consider its relation to known objects in 
both the Cygnus region and the rest of the Galaxy. Second, we find that a simple hadronic model can easily 
accommodate Milagro's flux measurement (which is at a single energy), as well as other existing observations 
spanning nearly seven orders of magnitude in gamma-ray energy. Third, since a hadronic gamma-ray spectrum 
necessitates an accompanying TeV neutrino flux, we show that IceCube observations may provide the first direct 
evidence of a Galactic cosmic-ray accelerator. 



I. INTRODUCTION 

The Cygnus region is one of the most prominent features 
of nearly every Galactic skymap, across many orders of mag- 
nitude in energy. In gamma rays, sources have been known 
to populate this region since the time of COS B fin. In 
subsequent years, EGRET discovered both diffuse 121 and 
point-like MeV-GeV emission [3], followed by the HEGRA 
observation of an unidentified TeV gamma-ray source [4], 
TeV J2032+4130. Recently, at even higher gamma-ray en- 
ergies, the Cygnus region has been discovered yet again by 
Milagro @], The gamma-ray flux measured by Milagro in- 
cludes a large (~ 15° x 10°) diffuse region and a new, uniden- 
tified TeV source, MGRO J2019+37 (see Fig. Q}. To clar- 
ify the relative positioning of these objects, Fig. [2] shows the 
regions containing TeV J2032+4130 and MGRO J20 19+37 
(which are separated by ~ 5°), as well as sources from the 
Third EGRET [3] and GeV la] catalogs (which all remain 
unidentified), overlaid on the diffuse emission measured by 
EGRET HH. 

For most TeV sources, particularly those unidentified at 
other wavelengths, it remains a mystery as to whether the ob- 
served gamma-ray flux is produced by hadronic or leptonic 
processes. Fortunately, neutrino telescopes offer hope for dis- 
criminating between these scenarios. If the gamma rays arise 
from the the decay of neutral pions (n° — > yy) produced in 
proton-proton scattering, it is well-established that a corre- 
sponding flux of neutrinos (produced in jt* decays) must also 
be present In contrast, no neutrinos result from the inverse 
Compton scattering (e~y — > ye~) of energetic electrons on 
ambient photons. 

HESS [9] has contributed significantly to TeV astrophysics 
by dramatically increasing the number of known southern-sky 
TeV sources (of interest HlOn to a km 3 Mediterranean neu- 
trino telescope filll ). Milagro's ability to survey the entire 
northern TeV sky is of great utility in making predictions 
for IceCube fl2tl . which, being located at the South Pole, is 
well-situated to observe upgoing muons initiated by neutrinos 
from northern-sky sources. The Milagro skymap [5] reveals a 
number of sources that will also be of interest to the ongoing 
and upcoming IceCube, MAGIC HI, VERITAS H, and 



GLAST lfl5ll projects, whose improved sensitivity over the 
previous generation of experiments will greatly enhance our 
understanding of Galactic cosmic -ray, gamma-ray, and neu- 
trino production. In particular, VERITAS is already surveying 
the Cygnus region in TeV gamma rays with excellent angular 
resolution and flux sensitivity. 

It is of advantage to understand this region in general and 
MGRO J20 19+37 in particular. Towards this goal, we will 
first examine the vicinity of this source to note possible identi- 
fied counterparts, which, among other things, give us concrete 
distances. Furthermore, in addition to Milagro, this region 
has been observed by other experiments in the past. We make 
use of these observations, particularly those of EGRET and 
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FIG. 1: The field near MGRO J2019+37. Shown are 3EG Cata- 
log sources, a GeV Catalog source, and potential counterparts to the 
gamma-ray sources. These sources are overlaid upon diffuse GeV 
emission observed by EGRET (white most intense). For details, 
please see the caption of Fig. [2] 
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FIG. 2: Gamma-ray sources and diffuse GeV emission in the Cygnus region. Shown are the sources discovered by Milagro (MGRO J2019+37) 
and HEGRA (TeV J2032+4130), along with their approximate (lc) error circles. The fitted extent of the Milagro source is comparable to the 
circle shown. Also shown are nearby Third EGRET (3EG) (compiled from > 100 MeV gamma rays) and GeV (> 1 GeV gamma rays) catalog 
sources (all at 95% confidence); as well as gamma-ray source candidates (points), the Cyg OB2 core (dashed circle), and the region of Fig. Q] 
(box). EGRET 4—10 GeV (point-source subtracted) diffuse emission (smoothed and scaled linearly from ~ 1 — 10 x 10~ 6 cm~ 2 s _I sr _1 , 
with white most intense) is also displayed. 



CASA-MIA, to construct a gamma-ray spectrum for MGRO 
J2019+37. Finally, we discuss the possibility of using Ice- 
Cube to determine whether the observed gamma-rays are the 
product of hadronic processes, which would indicate a site of 
cosmic-ray production. 



II. A TEV SOURCE OF UNKNOWN NATURE 

Most conspicuous in the Milagro view of the Cygnus re- 
gion is the new TeV source MGRO J2019+37 (see Fig. 3 of 
Ref. JH). While this source has been observed with very high 
significance (~ lie) at very high energies (median energy of 
12 TeV), no obvious multi- wavelength counterpart seems to 
be present. Analysis by Milagro does suggest that the TeV 
emission originates from either a single extended source or a 
combination of several unresolved point sources, fit by a 2D 
Gaussian of width = 0.32°. This implies a radius of 

'-(o§0(£)" 

which depends upon the unknown source distance (scaled to 
1 kpc for convenience). We present in Fig. Q] objects in the 
vicinity of MGRO J20 19+37 that merit further study. These 
objects (typically located at distances of ~ 1 — 4 kpc) can also 
be seen in Fig. [2] in the context of the greater Cygnus region, 



which (for the coordinate range displayed) is also nearly com- 
pletely covered by diffuse TeV emission |5|]. 

Of most immediate interest are two unidentified EGRET 
sources, 3EGJ20 16+3657 and 3EG J2021+3716 (constructed 
from gamma rays with energy > 100 MeV), and the GeV cat- 
alog source, GeV J2020+3658 (> 1 GeV gamma rays @]). 
These sources have received a fair amount of attention over 
the years in searches for potential counterparts. For 3EG 
J20 16+3657, a probable association with a blazar (of un- 
known redshift) was found 11611 . However, pair produc- 
tion on the extragalactic infrared background makes this an 
unlikely 12 TeV source [5]. 3EG J2021+3716 and GeV 
J2020+3658 have intriguing possible correlations with sev- 
eral Galactic objects. One of these is the pulsar wind nebula, 
PWN 75.2+0.1 [17], which, considering the growing number 
of known TeV PWNe (e.g., Ref. II 1811 ). is also one of the better 
candidates for MGRO J2019+37. 

In light of the recent possible discovery of a TeV source 
(HESS J1023-575) |QJ] coincident with the Wolf-Rayet 
eclipsing binary WR 20a [ 20], it is worthwhile to consider 
similar systems near MGRO J2019+37, some of which have 
been examined as possible counterparts to the EGRET sources 
(e.g., Ref.j2l). Shown in Fig. Q] are the Wolf-Rayet stars 
WR 137 irHTwR 138 |H, WR 141 fH, WR 142 fH, 
and HD 228766 [24]. Also displayed is V382 Cyg, a mas- 
sive eclipsing binary with a significant rate of mass loss l25ll . 
In fact, many of these systems are known to have large mass 
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loss rates (> 10~ 5 M Q yr _1 j26ll ). which may power shocks 
that accelerate cosmic rays. The reader is encouraged to see 
Ref. Il27ll (and references therein) for details of gamma-ray 
modeling in such environments. In a similar vein, the open 
star cluster Berkeley 87 has been proposed as a source of 
gamma rays J28j]. 

With the class of TeV supernova remnants already having 
several prominent members II29I1 . such an object in Cygnus 
would be a compelling candidate. However, a search of the 
area immediately around MGRO J20 19+37 yields only SNR 
074.8+00.6 (likely just an HII region (U) and the quite dis- 
tant (12 kpc (H) SNR 074.9+01.2 (an unlikely association). 
The high energy source closest to the quoted Milagro position, 
IGR J20188+3647, is seen only in the 17 - 30 keV band 03- 
Also worth noting is Cygnus Rift, a molecular cloud com- 
plex a few kpc in extent which runs east-west through this 
region [31]. 



cated near this area, perhaps the foremost being the micro- 
quasar SS 433 and its SNR, W50, from which no TeV emis- 
sion was discovered by HEGRA ll4lll . Other systems worth 
noting are SNR G40.5-0.5 J4I, HMXB 4U 1909+07 H, 
candidate SNR 041.3-01.3 114411 . and a prominent molecular 
cloud complex j45ll . Clearly, more observations are needed in 
order to make any further conclusions; however, this appears 
to be another exciting region for future studies in high-energy 
astrophysics. 

Since every source observed by Milagro, except for 
Mrk 421, is located within the Galactic Plane, these observa- 
tions strongly suggest that MGRO J20 19+37 is itself situated 
within the Milky Way. We will proceed with further analysis 
following this assumption. 



IV. GAMMA RAYS AND COSMIC RAYS 



III. WHAT ELSE IS OUT THERE? 

Thus far we have only entertained Galactic objects as po- 
tential counterparts to MGRO J20 19+37. When we exam- 
ine the locations of other sources in the Milagro map of the 
northern sky [5], additional support is found for this scenario. 
In addition to MGRO J20 19+37 and the large region of dif- 
fuse emission from Cygnus, Milagro has previously identified 
the Crab nebula and the blazar Markarian 421 13311 . two clas- 
sic TeV sources. Milagro has also confirmed TeV emission 
from the vicinity of the well-studied 10, HU, but poorly under- 
stood, TeV J2032+4130, which appears distinct from MGRO 
J2019+37. Assuming that these sources are ~ few kpc from 
Earth, their angular separation implies that the distance be- 
tween them is at least a few hundred pc, which suggests that 
they are unrelated. 

Since Milagro has, essentially, the same field of view as Ice- 
Cube, it is also useful to note any other objects which may be 
of interest to neutrino (and gamma-ray) astronomy. We con- 
sider the brightest regions in the Milagro skymap, for which 
identifications seem to be evident, but which were not noted in 
Ref. |5j|. Of course, since these regions have not yet been ad- 
dressed in particular, further analysis by the Milagro Collabo- 
ration will provide the final word. Besides the Cygnus region, 
along the Galactic Plane are several other regions of possi- 
ble diffuse emission, which (if confirmed) would help us to 
further understand the "TeV excess" IB5I1 problem that arose 
from the initial Milagro discovery of Galactic Plane emis- 
sion OHim and remains an outstanding concern [5]. These 
are located around Galactic longitudes I rj 50° and I « 35° 
and roughly coincide with the tangents of the Milky Way's 
Sagittarius and Scutum arms JH HH, respectively. Confir- 
mation of such emission would have tremendous value in the 
study of cosmic rays and their interaction with matter in the 
Galaxy @- 

Also in the Plane, near I — 40°, is a point of > 6(7 signifi- 
cance nearly directly coincident with GeV J 1907+0557 ]6[]. It 
is intriguing that a Whipple observation of this region showed 
an excess of gamma-ray-like events over the entire field of 
view, which was deemed unlikely to be due to a bright, ex- 
tended source 114011 . A number of interesting objects are lo- 



Additional insight can be obtained by constructing a 
gamma-ray spectrum based on the Milagro measurement. Mi- 
lagro reported a gamma-ray flux from MGRO J20 19+37 of 
E 2 d<P/dE = (3.49±0.47 sta t±1.05 sys ) x 10" 12 TeV cm^s" 1 
at a median energy of 12 TeV (see Fig. [3]), with an undeter- 
mined (but subdominant) contribution from the surrounding 
diffuse emission. This detection, at such a high energy, is 
quite useful in creating a spectrum, however, we also need 
additional data at higher and lower energies. The two nearby 
(unidentified) EGRET sources may provide information at the 
low end. Both have similar spectra in the MeV-GeV range 0]; 
however, as 3EG J2021+3716 is more likely [46] to be asso- 
ciated with GeV J2020+3658, we will consider its possible 
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FIG. 3: Data and possible hadronic spectra for MGRO J2019+37. 
Shown are the Milagro measurement at 12 TeV (diamond), the 
EGRET spectrum for 3EG J2021+3716 (circles), the upper limit 
from Whipple (0.3 TeV), and our inferred upper limit from CASA- 
MIA (115 TeV). Also shown are hadronic fits to the data, assuming 
£/7 2 ' 35 ( u PP er ) an d Ep 2 (lower) source proton spectra. The region 
above 1 TeV is of greatest interest to neutrino astronomy. 
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association with the Milagro source. Near 1 TeV, Whipple 
observations place upper limits on emission from the region 
near 3EG J2021+3716 (see Fig.[3]i, with a stronger limit from 
the vicinity of PWN 75.2+0.1 114011 (not shown). Emission 
from Berkeley 87 is also constrained from HEGRA observa- 
tions [47] (also not shown). 

At energies higher than the Milagro measurement, obser- 
vational limits are crucial in determining the spectral cutoff, 
which is of great importance to cosmic-ray and neutrino stud- 
ies IIOII . Fortunately, CASA-MIA observations can provide 
valuable limits. While a previous all-sky limit is not very con- 
straining j48l], later searches had greatly enhanced statistics 
allowing for an improvement of more than an order of mag- 
nitude [14911 . Although there is yet no published limit specif- 
ically for the MGRO J2019+37 region, we can infer a limit 
based upon observations of the region containing Cygnus X- 
3 14^1. as there were no significant sources in its general vicin- 
itv Il49l [50h . The integral flux limit from CASA-MIA at the 
position of Cygnus X-3 is <&(£ > 115 TeV) < 6.3 x 10~ 15 
photons cm~ 2 s , given above the median detector energy 
in order to reduce dependence on the assumed spectral in- 
dex 14911 . We infer a similar limit for the region near MGRO 
J20 19+37 by conservatively rounding the Cygnus X-3 result 
up to 10 _14 cm _2 s _1 and by treating the integral as Ed'P/dE, 
giving an upper limit of E 2 d<P/dE = lQ- n TeV cm" 2 s" 1 at 
1 15 TeV for the MGRO J2019+37 region (see Fig. E), which 
(even if somewhat higher) gives meaningful constraints. 

As the properties of this source (e.g., magnetic fields) that 
would determine a leptonic spectrum are quite unknown, we 
will instead construct a simple, entirely hadronic spectrum 
(which requires fewer assumptions), not attempting to include 
contributions from primary or secondary electrons. We use the 
parametrization for gamma rays resulting from p-p scatter- 
ing of Kelner, Aharonian and Bugayov 13 lfl - The gamma-ray 
spectrum resulting from a source distribution of protons of the 
form d<t> p /dE p =A p E p a exp [- (E p /Ep Ut )] , interacting with 
ambient protons of density tin , is 

d® 7 r , ^ d<Z>„ (E y \ dE„ 

iir cnH i a ^ (Ep) ^A^ Ep )^' (2) 

where Oi ne i(Ep) is the inelastic p-p cross section, c is the 
speed of light, and the function Fy(E 7 /E p , E p ) (given by 
Eq. (58) of Ref. [51]) determines the number of photons pro- 
duced per energy interval per scattering. We use Eq. (01 for 
Ey > 1 TeV and the c) -function approximation of Ref. 115 ill at 
lower energies to produce our gamma-ray spectra. 

We first assume that there is an association between 3EG 
J2021+3716 and MGRO J2019+37. Using an £~ 2 - 35 proton 
spectrum (with £™' = 1000 TeV) and normalizing to the Mi- 
lagro measurement at 12 TeV, we find the upper spectrum of 
Fig. |3] This spectrum gives a reasonable (considering the un- 
certainties involved) fit to the EGRET data and, importantly, 
satisfies the CASA-MIA limit. While the distance (^) to the 
source and ambient proton density («#) are both unknown, 
normalizing to the Milagro flux allows us to find the total nec- 
essary energy injected into cosmic -ray protons as 



which, for comparison, is similar to the total explosion energy 
of a typical core-collapse supernova [52] when @ ~ few kpc. 

It is quite possible that neither of the EGRET sources are 
at all related to MGRO J2019+37. In any case, EGRET (as 
well as Whipple and HEGRA) constrains the spectrum from 
extrapolating to lower energies too steeply, if only for the sim- 
ple reason that another source would have been seen. We 
also consider an Ej 2 proton spectrum (with Ep Ut — 500 TeV), 
which gives the lower spectrum in Fig. [3] While safely be- 
low the low-energy observations, the CASA-MIA limit again 
forces a cutoff at high energies. This scenario requires ~ 
10 times less input cosmic -ray energy, due to the lower re- 
quired GeV gamma-ray flux. 

The most economical approach towards explaining the Mi- 
lagro observation would be to have the gamma-ray spectrum 
peak (in our E 2 plot) at ~ 10 TeV and then sharply fall off (to 
satisfy CASA-MIA). Such a spectrum has been observed from 
the PWN Vela X by HESS [Dl, which may have a hadronic 
origin [53]. If PWN 75.2+0.1 has a spectrum of this form, it 
would easily satisfy the stronger Whipple limit. However, at 
distance of ~ 10 kpc [17], it would need to be more than an 
order of magnitude more energetic than Vela X to account for 
the Milagro measurement. 



V. CYGNUS AND ICECUBE 

While a spectral analysis gives important guidance, the 
most direct way to discern the nature of a TeV source is 
through neutrino observations. In high energy p-p scatte ring , 
7t + , 7t~, and 7T° are produced in roughly equal numbers f54fl . 
Gamma rays result from the decay n° — > yy, while neutrinos 
arise from the 7t + — > ^i+Vu — > e + v v y e v p _ and n~ — » \i ~ v M — » 
e^VfiVeVfi decay channels. The ratio of neutrinos to photons 
is found by considering the initial neutrino flavor ratio from 
charged pion decay, v e : : V T = 1 : 2 : 0, which is trans- 
formed into 1 : 1 : 1 by vacuum neutrino oscillations before 
reaching Earth. In neutrino telescopes, neutrinos and antineu- 
trinos are practically indistinguishable, so we typically con- 
sider the sum, v + V. Each photon from K° decay then corre- 
sponds to one neutrino of each flavor (Ny = N Ve = N VlA = N Vt ). 
The typical neutrino energy is ~ 1/2 of the gamma-ray energy 
from 7T° decay. The v + V spectrum is thus shifted, relative to a 
power-law gamma-ray spectrum of form d<Py/dEy = <j>yEy T , 
as 

^ = (l/2) r - 1 «M? v - r = </> v £ v - r , (4) 

where each neutrino flavor is treated separately. Our interest 
will be limited to the + flux, which can produce ultra- 
relativistic muons in charged-current interactions within the 
Antarctic ice cap. 

Knowledge of the location of gamma-ray sources a priori 
allows for neutrino searches without having to pay the trials 
factor associated with an undirected, all-sky search. Knowl- 
edge of source spectra then allows for neutrino flux predic- 
tions, which can be compared with observations. For MGRO 
J2019+37, we can use the spectra shown in Fig.[3]to derive the 
associated neutrino flux and the expected event rate of TeV 
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FIG. 4: Integrated (v^ + v f i)-induced muon rates from MGRO 
J20 19+37 above a given muon energy within IceCube for one year. 
These rates result from an £^ 2,35 source proton spectrum (see 
Fig. f3). The E~ 2 fit yields a nearly identical result. The shaded 
region shows the expected atmospheric background in a 3 deg 2 bin. 

muons in IceCube. Following Ref. 15111 . we parametrize each 
gamma-ray spectrum with a fit of the form 

§^A r £ f V(W<) 1/2 , (5 ) 
dE 7 

with Ay chosen to fit the gamma-ray spectrum at 1 TeV. For the 
E p 235 proton spectrum, we use j3 = 2.2 and Ey Ut = 45 TeV; 
while we find J3 = 1 .9 and Ey at = 20 TeV well fits the gamma 

rays resulting from the E~ 2 input spectrum (above 1 TeV). 
Using the methods detailed in Ref. [ 10], we then calculate the 
expected spectrum of (v^ + v^)-induced muons (as a func- 
tion of the observed muon energy entering or originating in 
the detector) for IceCube, as seen in Fig.|U The uncertainties 
in these results are at or below that of the Milagro flux mea- 
surement itself. This result is compared to the background 
arising from the atmospheric neutrino spectrum (which falls 
off more steeply with energy) in a 3 deg 2 bin, which roughly 
corresponds to the angular resolution of IceCube. Note that 
we only consider muons with E^ > 1 TeV, which are more 
likely to have an accurately measured energy and direction 
than less energetic muons. This is important for effectively 
discriminating between the expected signal and background. 

With such a similar gamma-ray spectrum above 1 TeV, the 
Ep 2 result is practically indistinguishable from the Ep 2 " 35 
case displayed. In fact, since the Milagro measurement is at 
such a high energy, the expected muon rate is only weakly de- 
pendent on the spectral index, with the location of the cutoff 
being somewhat more important. This rate satisfies the con- 
straints previously derived in Ref. II 1011 for the Cygnus region 
and is consistent with the results of Ref. [55], which we rec- 
ommend the reader to examine for an alternative approach. 



We find that IceCube can detect neutrinos from the MGRO 
J20 19+37 source, if indeed it is a cosmic-ray accelerator. It 
is worth emphasizing that not a single high-energy neutrino 
source has yet been detected. While the number of events 
will not be large, a signal could be separated from background 
within several years, especially taking advantage of the rapid 
increase in signal to background probability with increasing 
measured muon energy. 

VI. PROSPECTS AND CONCLUSIONS 

Milagro has discovered a unique and very interesting ob- 
ject, MGRO J2019+37, in the Cygnus region of the Galaxy: a 
source bright in ~ 12 TeV gamma rays, apparently extended 
in size, and without clear counterparts in other wavelengths. 
At present, both the nature and distance of this source remain 
unknown. A near-term experimental approach that utilizes 
a variety of techniques is immediately suggested. GLAST 
observations of MeV-GeV gamma rays should determine the 
validity of the EGRET discovery of several distinct sources 
in this region and eliminate any source confusion that may 
be present (as evidenced by the location of the GeV catalog 
source). Whatever objects are actually present, their MeV- 
Ge V spectra will be measured with greater accuracy. This will 
be important in discerning whether a connection exists with 
the Milagro measurement at 12 TeV (e.g., the upper spectrum 
of Fig. 

Observations in the GeV- TeV range by VERITAS and 
MAGIC will provide a clear picture of the morphology of 
the emission and conclusive evidence as to whether single or 
multiple TeV sources are present. Measurements of the shape 
(steep or shallow) of the gamma-ray spectrum will allow for 
improved source modeling and calculations of the associated 
neutrino flux. Experiments like those of the Milagro and Ti- 
bet 115 611 groups will play a significant role by prospecting for 
sources at the highest energies, and by mapping out diffuse 
emission. Sometimes lost in discussions of whether hadronic 
TeV sources exist is the fact that the cosmic -ray protons ob- 
served in abundance at Earth must be produced somewhere in 
the Galaxy. Given the estimates of source and interstellar gas 
densities, pion-producing collisions should be expected. 

To probe the nature of this object, we have pursued three 
lines of attack. First, we considered in detail the relation of 
this object to others in the Cygnus region, as well as the rest 
of the Galaxy. The Cygnus region is a site of intense, ongoing 
star formation that is dense in observed and candidate gamma- 
ray sources, as well as diffuse emission. The source distance 
is still unknown, but could be ~ 1 — 2 kpc if associated with 
other objects in this region lf57ll . Second, while the Milagro 
measurement is only in a single energy bin, we showed how 
the entire gamma-ray spectrum is constrained with existing 
observations. We find that a simple hadronic model can eas- 
ily fit EGRET data at lower energies, as well as the inferred 
CASA-MIA limit at higher energies. Third, the hypothesis of 
a hadronic spectrum can be decisively tested by the accompa- 
nying neutrino flux, and we showed that IceCube should be 
able to observe this object in several years of operation. The 
detection of muons initiated by TeV neutrinos would thus au- 
thenticate a Galactic cosmic-ray accelerator. 
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